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f*^ , In this talk, I review the important role played by plasma instabilities in the ther- 

f— *^ ■ malization of quark-gluon plasmas at very high energy. [Conference talk presented 

' at Strong and Electroweak Mattter 2004, Helsinki, Finland, June 16-19.] 

o ' 
o 



1. Introduction 



Here is a basic question: What is the (local) thermalization time for quark- 
gluon plasmas (QGPs) in heavy ion collisions? That's a difficult question, 
so let's ask a simpler one: What is it for arbitrarily high energy collisions, 
, where the running coupling constant is small, a s <C 1? Even that turns 

out to be complicated, so let me focus on an even simpler version, first 
posed by Baier, Mueller, Schiff and Son: 1 How does that time depend on 
a s ? In particular, in the saturation picture of heavy-ion collisions, what is 



' the exponent n in the relation 

a, 



Before discussing this question in more detail, I wish to make a general 
side comment about plasma physics. Plasma physics is complicated! This 
is made abundantly clear simply by looking at pictures of various plasma 
phenomena, such as the image in Fig. 1 of a solar coronal filament from 
NASA's TRACE satellite. Theoretical discussions of quark-gluon plasmas, 
however, are generally much less complicated. There arc several reasons 
why such discussions can usually avoid the full complication of traditional 
plasma physics. First, much of the theoretical literature discusses QGPs 
that are at or very close to thermal equilibrium. The physics of plasmas 
near global thermal equilibrium is much less complicated than the physics of 
non-equilibrium situations. Because electromagnetic interactions are long- 
ranged, traditional electromagnetic plasmas can be very complicated even 



1 
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when they are in local thermal equilibrium. Macroscopic currents in one 
region of the plasma can interact magnetically with other currents in other 
regions, over tremendous distance scales, creating complicated structures 
like Fig. 1. Non-Abelian plasmas, however, are somewhat different. From 
theoretical studies of the equilibrium properties of such plasmas, we know 
that the non-Abelian interactions cause magnetic confinement over dis- 
tances of order l/(g 2 T). It is reasonable to assume that, even dynamically, 
color magnetic fields cannot exists on distance scales larger than the con- 
finement length. So, unlike traditional electromagnetic plasmas, there are 
no large- distance magnetic fields. As far as the color degrees of freedom 
are concerned, the long-distance effective theory of a non-Abelian plasma 
is hydrodynamics rather than magneto-hydrodynamics. 




Figure 1. Image of a solar coronal filament from NASA's TRACE satellite, from 
(http:/ /antwrp. gsfc.nasa.gov/apod/ap000809. html). 

The color magnetic fields can play a role on small scales. But the full 
complications of plasma physics might be ignored on small distance scales if 
the relevant physics on those scales is weakly interacting. This was the pro- 
posal of the original bottom-up scenario for thermalization of quark-gluon 
plasmas* However, as we shall discuss, even at small distance scales, there 
can be plasma instabilities. These instabilities cause the growth of non- 
perturbatively large magnetic fields, bringing in some of the complicated 
non-linear physics of traditional plasma physics. 
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In the original bottom- up thermalization scenario of Baier et ai, the 
starting assumption is the saturation scenario, which assumes that the for- 
mation of the quark-gluon plasma starts at times t ~ Q" 1 dominated by 
gluons with momenta p ~ Q s , where the scale Q B is known as the saturation 
scale. In this talk, these gluons will be referred to as "hard" gluons, since 
we will soon be discussing even softer momentum scales. This situation 
is depicted in Fig. la. The occupation numbers of each gluon mode are 
initially non-perturbative, with f(x,p) ~ l/a s , where / is the phase-space 
density. As the system expands 1-dimensionally immediately after the colli- 
sion, the density per unit volume decreases, and one might therefore expect 
the hard gluon interactions to become more perturbative. To understand 
what happens next, let's ignore these perturbative interactions for the mo- 
ment and think about free expansion. As the nuclear pancakes fly apart 
after the collision, the gluons, which started in the center, will separate 
themselves in space according to the z components of their velocities, as 
Az ~ v z t. This is depicted in Fig. lb. The gluons left in the central region 
will be those with small v g , as shown in Fig. lc. As a result, the momentum 
distribution in that central region will have an anisotropic pancake shape, 
as shown in Fig. Id. In other regions, the momentum distribution is simply 
a boosted version of this — that is, it looks the same if one works in the local 
rest frame. 
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Figure 2. Approximately free expansion at relatively early stages of the bottom-up 
scenario. 



It should be emphasized that this anisotropic momentum distribution 
has nothing to do with the usual elliptic flow distribution measured in non- 
central heavy ion collisions. Fig. Id is a statement of the local momentum 
distribution in the local fluid rest frame at very early times, before collisions 
have brought the system to local equilibrium. Elliptic flow is instead a 
measure of the net fluid flow on large scales of a system that has quite 
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possibly come to local thermal equilibrium (and so locally has isotropic 
momentum distributions in the local fluid rest frame). 

In the original bottom-up scenario, equilibration of the plasma was 
assumed to occur through individual 2-body collisions between particles 
(with some LPM effect thrown in). In the first stage of the scenario, 
1 Q s t <C a~ 5 / 2 , the important processes were small angle scatter- 
ing, which slightly widens the hard particle distribution in p z , and soft 
Bremsstrahlung from colliding hard particles, which creates soft gluons 
with momenta k <C Q s . In the second stage, a~ 5 / 2 <C Q s t <C a -13 / 5 , 
these soft gluons come to dominate the number density of particles, but 
the hard gluons still dominate the energy density. Collisions between the 
soft gluons cause the soft gluons to thermalize. Finally, the hard particles 
begin to lose energy by Bremsstrahlung plus cascading, as shown in Fig. 1. 
In this scenario, local thermalization became complete at Q s t ~ a -13 / 5 . 



Figure 3. Cascading process for energy transfer from hard particles to soft sector in the 
original bottom-up scenario. Straight lines denote hard gluons, p ~ Q B . 

The original bottom-up scenario overlooked the possibility that collec- 
tive processes (as opposed to sequences of individual collisions) could play 
an important role in the equilibration of the plasma. In the case at hand, 
these collective processes are related to the appearance of plasma instabil- 
ities in the analysis of the equilibration of the quark-gluon plasma. 

2. Plasma Instabilities 

The hero of this story is Stan Mrowczyhski 2 ' 3,4 who over the years has been 
the major proponent of the idea that plasma instabilities are important for 
the equilibration of the quark-gluon plasma. The application of this idea 
to bottom-up thermalization was made by myself, Jonathan Lenaghan and 
Guy Moore? A selection of other folks who have considered the idea past 
and present include Refs. 6-9 and the work by Romatschke and Strickland* 
which was reported on at this conference. 
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Let me start with a slightly formal explanation of the origin of plasma 
instabilities, and I will give a more physical picture afterward. Imagine 
calculating the self-energy H(lu, k) for a particle moving through the plasma. 
The self-energy represents the effect on the particle of forward-scattering off 
of other particles in the medium, as in Fig. 2a, which one can alternatively 
calculate as in Fig. 2b using one's favorite formalism for field theory in a 
medium, or more simply calculate using linearized kinetic theory. (Here, the 
straight lines represent hard particles, which are gluons in the bottom-up 
scenario.) The result is the same either way. Generically, if the momentum 
distribution f(p) of hard particles is anisotropic, 21 one finds that there are 
negative eigenvalues of IP J (0, k). (See, for example, the general arguments 
in Ref. 5.) Such negative eigenvalues indicate instabilities at small k, which 
means exponentially growing soft gauge fields. 




(a) (b) 

Figure 4. Self-energy of soft modes due to hard particles. 

Let me give an analogy from scalar </> 4 theory at finite temperature. 
Imagine integrating out the hard particles to get an effective thermal self- 
energy II of order AT 2 . The effective linearized equation of motion is then 

uj 2 ^^ (k 2 + m 2 + n)cj) = (k 2 + m 2 cfi )(t). (2) 

Let's set mo = to improve the analogy, since gluons do not have any 
intrinsic mass. If we were in a situation where k 2 + H were less than 
zero, then there would be solutions with u> = ±iT pure imaginary, which 
would lead to exponentially growing solutions to the linearized equation. 
Alternatively, think about the effective potential, 

V eS (0) = m 2 ff </> 2 + A0 4 = U<j) 2 + A0 4 . (3) 



a Here and throughout, I always assume / is parity symmetric: /(— p) = /(p). 
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If II < 0, as happens in some multi-scalar theories} 1 then the potential 
looks like a double-well potential: the naive vacuum <p = is unstable; 
there is exponential growth from <fi = of modes with k < (—II) 1 / 2 ; and 
the growth stops once <j> becomes non-perturbatively large. 

Now let's turn to a physical picture of the instability, which I will adapt 
from Refs. 2,5. For simplicity, imagine two inter-penetrating, homogeneous 
streams of charged hard particles, one going up the page and one going 
down, which I'll call the ±x directions. Now also imagine that, due to 
some fluctuation, there is a very tiny seed magnetic field of the form B = 
Be y cos(fcz), as shown in Fig. 2a. Here, crosses denote magnetic fields 
pointing into the page, and dots fields pointing out of the page. Using the 
right-hand rule, you can check that the magnetic fields bend the trajectories 
of positively charged particles in the directions shown. This then focuses 
the net downward and upward currents into different channels, as shown in 
Fig. 2b. Again using the right-hand rule, one finds that these currents in 
turn create magnetic fields that add to the original seed field. With bigger 
fields, the effect becomes more pronounced, and the fields continue to grow 
through this mechanism. This is the Weibel instability. 




We get a seemingly contradictory picture of what happens if we consid- 
ers hard particles which move in other directions, such as in Fig. 2. Fol- 
lowing how these particles are by the seed magnetic field, we find they are 
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directed a little more upward in some regions and a little more downward 
in others. This results in a net current as shown in the picture, but these 
currents create magnetic fields which oppose the seed field. Let p be the 
original momentum of the hard particle and k the wave number of the soft 
magnetic field fluctuation. What happens is that particles with p ■ k ~ 
get trapped in valleys, as shown in Fig. 2b, and give the de-stabilizing effect 
discussed earlier. Other particles, with p k 0, are "untrapped" as in Fig. 
2a, and give a stabilizing effect. For isotropic f(p), these two contributions 
turn out to cancel, giving n y (0, fc) = 0. 




Now, instead of an isotropic hard particle distribution f(p), think of 
the one depicted earlier in Fig. Id. For k in the z direction, we will get 
a relatively smaller percentage of particles with p ■ k ~ than we would 
with an isotropic distribution, and so the net effect will be stabilizing. On 
the other hand, for k in the _L direction, we will get a relatively larger 
percentage of particles with p ■ k ~ than we would with an isotropic 
distribution, and so the net effect will be de-stabilizing. For k at significant 
angles to the axis, most of the particles with have p ■ k jk 0, and we will 
tend to get stability. The moral of this story, to which we will later return, 
is that the unstable modes associated with this distribution have k pointing 
very close to the z axis. 



February 2, 2008 0:46 Proceedings Trim Size: 9in x 6in 



sewm04 



8 

3. When does the growth stop? 

As in the scalar analogy discussed earlier, the growth of instabilities should 
stop when the soft fields become non-perturbatively large. From considering 
covariant derivatives D = d — igA, the effects of soft fields A become non- 
perturbative when A ~ d/g. There are two possibilities for the momentum 
scale associated with the derivatice d. The first, a possibility for both 
QCD and QED, is that growth stops when the effects of the soft fields 
on hard particles becomes non-perturbative, which will happen when A ~ 
Phard/ff and corresponds to the trajectories of hard particles being bent 
dramatically from straight lines. The second possibility, which cannot occur 
in QED, is that growth stops when the non-Abelian self- interactions of 
the soft fields become non-perturbative. This corresponds to A ~ k so f t /g. 
Whenever there is a significant separation k -C p between soft and hard 
physics (as in the bottom-up scenario), these two possibilities correspond 
to parametrically different scales for A. 

Jonathan Lcnaghan and I conjecture 12 that growing QCD instabilities 
"abelianize." That is, the growth stops when A ~ k so f t /g ~ |n| 1//2 /g, 
just as in QED. That in turn suggests that the complicated stuff that 
happens afterward is closely related to the mainstream plasma physics of 
(collisionless) relativistic QED plasmas. What follows is a summary of 
suggestive arguments that we make for abelianization. 

Start with the general HTL effective action for anisotropic f(p), which 
I adapt from Mrowczyriski, Rebhan, and Strickland 13 Shcematically, 

S eS = - [ F 2 -g 2 f f(p)W 2 , (4) 

J x J xp 

W = W a (x,p) = ^F fia . (5) 

Now imagine finding the effective potential V e s by looking at S e s for time- 
independent configurations in ^4o=0 gauge. There is a problem, which is 
that (p ■ D)^ 1 is a complicated, non-local operator. But now recall that, 
for the hard particle distribution of Fig. Id, the typical unstable modes 
have k pointing very close the z direction. Inspired by this, let's ignore 
k±_ altogether and consider configurations A = A(z) depending on z only. 
There is then an amazing simplification, noted by Blaizot and Iancu: 14 W 
given by (5) then turns out to be linear in A. As a result, the HTL term 
g 2 J fW 2 in the effective action (4) is then quadratic in A. That means 
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that it consists of nothing but the HTL self-energy, so that 

S eS [A(z)] = - [ F 2 ~ [ AUA. (6) 

J x J X 

The IT here is for fc's proportional to z. It is also transverse, corresponding 
to II_i_(0, kz). This turns out not to depend on the magnitude k of k. I will 
define the constant — /i 2 = H±(0,kz). The effective potential from (6) is 
then 

V cS [A(z)] = B 2 -^ 2 A\, (7) 

where B is the full, non-Abelian magnetic field, which includes non-linear 
terms in A. Now consider the arbitrarily soft limit k — ► 0, where the 
derivative term in B vanishes and only the non-Abelian commutator term 
survives: 

V cS ^[A x ,A y ] 2 -^(Al + Al). (8) 

What does this potential look like? As an example, let's consider just two 
of the degrees of freedom: color 1 of A x and color 2 of A y . A plot of the 
resulting potential is shown in Fig. 3. As one moves away from A = 0, the 
potential bends down due to the —fjL 2 A 2 term. If both colors are present, it 
then later bends back up again due to the quartic [A x , A y } 2 term. However, 
if the configuration is single-colored (A^p = or Ay 2 ^ — 0), the commutator 
vanishes, and the potential continues to run away. If we were to roll a ball 
from the origin in such a potential, it would eventually roll away down into 
one of the Abelian directions [A x , A y ] = —0. 




Figure 7. A picture of the effective potential. 

The above argument suggests the following conjecture, which is the con- 
clusion of my talk: The growth of Wcibcl instabilities eventually "abelian- 
izes" the soft gauge fields into the maximal Abelian subgroup of the gauge 
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group. For SU(2) gauge theory, this would give traditional U(l) plasma 
physics. For SU(3) gauge theory, it would give U(l)xU(l) plasma physics, 
corresponding to two copies of Abelian electromagnetism. Further discus- 
sion, including corroboration from numerical simulation results of a related 
toy model, may be found in Ref. 12. 
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